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We speak of chaos in quantum systems if the statistical properties of the eigenvalue spectrum 
coincide with predictions of random-matrix theory. Chaos is a typical feature of atomic nuclei 
and other self-bound Fermi systems. How can the existence of chaos be reconciled with the known 
dynamical features of spherical nuclei? Such nuclei are described by the shell model (a mean-field 
theory) plus a residual interaction. We approach the question by using a statistical approach (the 
two-body random ensemble): The matrix elements of the residual interaction are taken to be 
random variables. We show that chaos is a generic feature of the ensemble and display some of its 
properties, emphasizing those which differ from standard random-matrix theory. In particular, 
we display the existence of correlations among spectra carrying different quantum numbers. These 
are subject to experimental verification. 



PACS numbers: 21.60.Cs,24.60.Lz,21.10.Hw,24.60.Ky 



I. RANDOM MATRICES AND CHAOS 



A. Introduction 



This paper contains a colloquium talk on the appli- 
cation of random matrices to nuclear spectroscopy. We 
discuss the origin of random-matrix theory (RMT), some 
of its predictions, and its relation to classical chaos. We 
show that the predictions of RMT often agree very well 
with spectroscopic data in nuclei. In the main part of the 
paper, we address the question: How can this success of 
RMT be reconciled with our knowledge of the dynami- 
cal behavior of nuclei as embodied in the nuclear shell 
model? Although actually formulated for nuclei, we be- 
lieve that our arguments apply likewise in modified form 
to other Fermi systems like atoms and molecules. 

Near the ground state, the spectra of self-bound Fermi 
systems are essentially discrete. As an example, Fig. [T] 
shows the measured spectrum of the nucleus 19 0. Up 
to the threshold for decay into a neutron (n) and 18 
(the "neutron threshold" ) at an excitation energy of 3.957 
MeV, the levels shown decay only by gamma emission; 
the ground state is unstable against beta decay. The 
long lifetimes of these decay modes render the levels 
virtually discrete. Above neutron threshold, the levels 
acquire finite particle decay widths (which typically are 
much larger than those for beta and gamma decay) and 
are also seen as resonances in cross sections (like in the 
one for the reaction n + 18 0). The levels in Fig. [1] carry 
quantum numbers which reflect symmetries of the nu- 
clear Hamiltonian: Total spin J corresponds to rotational 
invariance, total isospin T to proton-neutron symme- 
try, parity P to mirror reflection symmetry. In heavy 
nuclei, the Coulomb interaction between protons breaks 
the isospin quantum number. In its general appearance, 
Fig. [1] is representative not only for nuclei but also for 
atoms, molecules, and atomic clusters. 
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FIG. 1 The spectrum of 19 0. Taken from Ref. (0). 



B. Regular Motion 

The excitation energies and widths for beta, gamma, 
and particle decay of the levels found experimentally can 
often be successfully described in terms of simple inte- 
grable dynamical models. In analogy to classical me- 
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chanics, we then speak of regular motion. A striking ex- 
ample is furnished by the existence of rotational bands. 
These are characterized by a spin/parity J p sequence 
+ ,2 + ,4 + ,... and excitation energies proportional to 
J {J + 1). They correspond to a rotation of the entire 
(non-spherical) nucleus about some axis. Fig. [2] shows 
two such bands in the nucleus 174 Hf. Another example for 
regular motion is provided by the independent-particle 
model: nucleons move independently in the mean field of 
the nucleus. This model corresponds to the most elemen- 
tary version of the nuclear shell model and accounts for 
the existence of "magic numbers" for neutrons and pro- 
tons (i.e., numbers where a major shell is filled and where 
the nuclear binding energy attains maxima) as well as for 
spins and parities of ground states of nuclei differing by 
one unit in proton or neutron number from closed-shell 
nuclei. We discuss the nuclear shell model in more detail 
in the next Section. 
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FIG. 2 Rotational bands in 174 Hf. Taken from Ref. (@). 



The first evidence for chaotic motion in nuclei came 
from spectroscopic data on levels at neutron threshold, 
i.e., rather far from the ground state. In general, it is 
exceedingly difficult to unambiguously identify the posi- 
tions, spins, and parities of such levels. Neutron time- 
of-flight spectroscopy offers the unique opportunity to 
acquire precise spectroscopic information in a range of 
excitation energies which is otherwise virtually inaccessi- 
ble. The neutron threshold in medium-mass and heavy 
nuclei typically occurs at an excitation energy of 8 to 
10 MeV (and not at w 4 MeV as in 19 0, see Fig. [Q). 
Measuring the total cross section for slow neutrons ver- 
sus kinetic energy, one observes resonances. Each of the 
resonances corresponds to a nuclear state at an excita- 
tion energy of around 8 or 10 MeV. The ground states 
of even-even nuclei have spin/parity + . Slow neutrons 
carry angular momentum zero and have spin 1/2. For 
even-even target nuclei, the states seen, therefore, all 
have spin/parity l/2 + . We show here not the earliest 
but some of the best data. In the 1970s, the Columbia 
group around Rainwater measured time-of-flight spectra 
of slow neutrons scattered on a number of heavy nuclei. 
For each target nucleus, they observed a sequence of up 
to 200 l/2 + levels. (That number was limited by the res- 
olution of the time-of-flight spectrometer.) Fig. [3] shows 
their data for the total neutron scattering cross section 
on 238 U versus neutron energy. We see a number of very 
narrow resonances (widths typically smaller than 1 eV) 
with typical spacings of 10 eV. The energies of the associ- 
ated levels were determined with the help of an i?-matrix 
multi-level analysis (0). 




C. Non-regular Motion 

The emphasis in this paper is on chaotic motion in nu- 
clei. The previous Subsection serves as a reminder that 
there is strong evidence for regular motion in nuclei, es- 
pecially in the vicinity of the ground state. We now turn 
to the equally strong evidence for chaotic motion. We put 
the argument in the form of a semi-historical narrative 
which will occupy the remainder of this Subsection. 



FIG. 3 Total cross section for scattering of neutrons by 238 U 
versus neutron energy. Taken from Ref. (|4|). 

Already in the 1930s, similar data (of much inferior 
quality) had led N. Bohr to formulate his "compound 
nucleus hypothesis" (|5|). Bohr argued that the existence 
of narrowly spaced narrow resonances is incompatible 
with independent-particle motion in the nucleus. Indeed, 
simple estimates using a central potential yield single- 
particle 8-wave level spacings around 1 MeV (and not 
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10 eV as shown in Fig. [3]), and s-wave decay widths for 
neutron-instable states of around 100 keV (and not < 
1 eV as shown in Fig. [3]). Bohr argued that the exis- 
tence of such narrow resonances could not be understood 
without assuming strong interactions between the inci- 
dent neutron and the nucleons in the target. A wooden 
toy model (shown in Fig. [4]) was to demonstrate this idea 
of the compound nucleus. The billiard balls represent 
nucleons, the queue indicates the kinetic energy of the 
incident neutron, and the trough simulates the attrac- 
tive mean field. Since the publication of N. Bohr's 1937 
paper and until the discovery of the nuclear shell model 
in 1949, the idea of independent-particle motion in the 
nucleus was virtually inacceptable in the nuclear physics 
community. And as we shall see, understanding chaos in 
nuclei is almost synonymous to reconciling N. Bohr's idea 
of the compound nucleus with the nuclear shell model. 




FIG. 4 Wooden toy model simulating N. Bohr's compound 
nucleus. Taken from Ref. (|6|). 



D. Random Matrices 

When Wigner (7) introduced random matrices into 
physics, he did not refer explicitly to N. Bohr's 1937 
paper (5). Nevertheless we believe that his work was 
motivated by and can be seen as the mathematical for- 
mulation of the idea of the compound nucleus. It is note- 
worthy that the work started in 1951, two years after the 
discovery of the shell model. 

How can one deal with a situation where the con- 
stituents of a quantum system interact strongly? We 
suppose that we deal with a non-integrable system with- 
out remaining symmetries. By that we mean that extant 
symmetries of the Hamiltonian like total spin and par- 
ity have been taken into account; the Hamiltonian has 
block-diagonal form; the individual blocks are free of fur- 
ther symmetries; we focus attention on one such block to 
which we refer for simplicity as the Hamiltonian. Is it 
possible to make generic statements about the spectrum 
and eigenvalues of such a system? Following Wigner (|7|) , 
we consider the matrix representation H^ v of the Hamil- 
tonian in Hilbert space. Here /i, v = 1, . . . , N and N ^> 1 
(we let TV — > oo at the end of the calculation). Nuclei are 



time-reversal invariant. Therefore, we can choose a rep- 
resentation where H^ v is real and symmetric. There are 
no further symmetries beyond rotational symmetry and 
parity. Now comes the decisive and unusual step: Rather 
than considering the individual Hamiltonian of the actual 
physical system, we study an entire ensemble of Hamil- 
tonians all having the same symmetry. It turns out that 
while we cannot make generic statements about the indi- 
vidual Hamiltonian, such statements are possible about 
almost all members of the ensemble. 

The ensemble is defined in such a way that it combines 
generality with the symmetry of H: It consists of real 
and symmetric matrices. To avoid a preferred direction 
in Hilbert space, it is chosen to be invariant under those 
transformations in Hilbert space which preserve the sym- 
metry of H. These are the orthogonal transformations. 
The ensemble is defined in terms of a probability density 
in matrix space. That density has the form 

/ TV traced 2 ) \ yj 

Cex P ( X 2 — ) 11 ^ v ^ 

The last factor is the product of the differentials of the in- 
dependent matrix elements. That factor is orthogonally 
invariant, and so is the trace of H 2 in the exponent. The 
Gaussian factor is introduced because a cutoff is needed 
to render the integrals over matrix space convergent. The 
ensemble is characterized by a single parameter A, which 
has the dimension of energy. This parameter determines 
the mean level spacing of the ensemble, i.e., the mean 
value of the distance between two neighboring eigenval- 
ues. We will be interested in the fluctuations of the spac- 
ings of neighboring levels around that mean value, and 
in the correlations between spacings of different pairs of 
neighboring levels. These are characterized by certain lo- 
cal fluctuation measures which are introduced below. All 
measures depend on the level spacing expressed in units 
of the mean level spacing and are, therefore, indepen- 
dent of the parameter A and altogether parameter-free. 
The factor C is a normalization factor. The factor N 
in the exponent (with N the matrix dimension) guaran- 
tees that the spectrum has finite range. Expression (pQ) 
defines the Gaussian orthogonal ensemble (GOE) of ran- 
dom matrices. Dyson (8) showed that there are three 
canonical random-matrix ensembles defined by symme- 
try: the GOE, the Gaussian unitary ensemble (GUE) 
(for time-reversal non-invariant systems), and the Gaus- 
sian symplectic ensemble (GSE) (for systems with half- 
integer spin which are not rotationally invariant). GUE, 
GOE, and GSE carry the Dyson parameters (3 = 1,2, 
and 4, respectively. 

While the symmetry of the Hamiltonian and, thus, the 
form of the invariant measure in expression (pQ) is dictated 
by quantum mechanics, the choice of the Gaussian cutoff 
is arbitrary and dictated by convenience. It has been 
shown, however, that an entire class of different choices 
(like, for instance, taking the trace of the fourth power 
of H in the exponent) yields the same results for the 



4 



local fluctuation measures. This important property is 
referred to as universality of the ensemble. It guarantees 
that the local fluctuation measures predicted by the GOE 
are generic. 

Random-matrix theory (RMT) only predicts ensem- 
ble averages of observables (calculated by integrating the 
observable over the ensemble (pQ)). But experimentally, 
we always deal with a single Hamiltonian, not an en- 
semble. The applicability of RMT to a single system is 
guaranteed by the property of ergodicity: For almost all 
members of the ensemble, the ensemble average of an 
observable (as given by RMT) is equal to the running 
average of that observable taken over the spectrum of 
a single member. Ergodicity guarantees that the quan- 
titative and parameter-free predictions of RMT can be 
meaningfully compared with experimental data. 



E. Fluctuation Measures 

The two fluctuation measures most frequently em- 
ployed in analyzing data are the nearest-neighbor spac- 
ing (NNS) distribution and the Dyson-Mehta or A3 
statistic. Fig. [5] gives the NNS distribution Pp(s), i.e., 
the probability distribution of spacings, versus s (the ac- 
tual spacing of two neighboring eigenvalues in units of 
their mean spacing) for Dyson's three canonical ensem- 
bles. The distributions are characterized by level repul- 
sion at short distances (Pp(s) oc s@ for small s) and a 
Gaussian falloff for large s. We note that if the sys- 
tem were integrable, all eigenvalues would carry differ- 
ent quantum numbers, and level repulsion were absent. 
The A3 statistic measures correlations between eigen- 
value spacings. Let M(E) denote the total number of 
levels below energy E. Clearly J\f(E) is a step function 
which increases by one unit as E passes an eigenvalue. 
The A3 statistic measures how well Af(E) is on average 
approximated by a straight line. It is defined by 

A 3 (L) = ^min a , 6 -jr J ° + (AT(E')-aE'-bj dE^ . 

%) 

Minimization with respect to a and b determines the best 
straight line. The angular brackets denote an average 
over the initial energy Eq and over the ensemble. All 
energies are in units of the mean level spacing, and A^(L) 
is, therefore, parameter free. If the spacings were totally 
uncorrelated, then A3 would be linear in L. It would 
grow much less fast with L if the spacings were correlated 
in such a way that a large spacing is always followed by a 
small one and vice versa. The actual behavior of As(L) 
is shown as a solid line in Fig. [2 ^(L) grows essentially 
logarithmically with L. This property is referred to as 
the "stiffness" of RMT spectra. 

For our later discussion, it is important to note that 
a GOE spectrum does not carry any information con- 
tent beyond the mean level spacing. This is obvious 
since any realization of H is obtained by drawing ran- 
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FIG. 5 The nearest-neighbor spacing distribution Pp(s) ver- 
sus s, the actual spacing in units of the mean level spacing, for 
the three canonical ensembles with /3 = 1, 2, 4 (solid, dashed, 
and dotted lines). [3 — 2 corresponds to the GOE. Taken from 
Ref. d). 

dom numbers from the Gaussian distribution (pp). Thus, 
it would be useless to apply spectroscopic analyses to 
such a spectrum. Put differently, a GOE Hamiltonian 
has N(N + l)/2 independent matrix elements. Counting 
shows that we need to measure all TV eigenvalues and all 
N eigenfunctions to reconstruct H. This is very different 
for typical dynamical systems where the Hamiltonian is 
known except for a small number of parameters. 



F. Quantum Chaos 

Suppose the fluctuation measures of an experimen- 
tal spectrum agree with GOE predictions. What does 
such a result imply physically? The understanding of 
this question was much advanced by the study of few- 
degrees-of-freedom quantum systems which are chaotic 
in the classical limit. The development culminated in the 
Bohigas-Giannoni-Schmit conjecture (10). It says: The 
spectral fluctuation properties of a quantum system which 
is chaotic in the classical limit, coincide with those of 
the canonical random-matrix ensemble that has the same 
symmetry. The authors of Ref. () 1 Qh supported their con- 
jecture by calculating numerically the NNS distribution 
and the A3 statistic for the Sinai billiard, a system which 
is time-reversal invariant and fully chaotic in the classical 
limit, and comparing the result with the GOE prediction. 
This is shown in Figs. [6] and [71 

Agreement like the one shown in Figs. [6] and has since 
been found in many other chaotic systems. Moreover, an 
analytical proof of the Bohigas-Giannoni-Schmit conjee- 



5 




0.5 1.0 1.5 2.0 x 2.5 



FIG. 6 The NNS distribution for the Sinai billiard (his- 
togram) and the GOE prediction (solid line). Taken from 

Ref. dla). 
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FIG. 7 The A3 for the Sinai billiard (histogram) and the 
GOE prediction (solid line). Taken from Ref. (flQh . 

ture has recently been published (11). That proof uses 
the semiclassical approximation and generic properties of 
chaotic trajectories. We are, thus, led to consider agree- 
ment between the spectral fluctuations of a physical sys- 
tem and RMT predictions as a signal for chaotic motion. 
It is in that sense that we will use the word chaos also 
in nuclei. This is done with a caveat: The semiclassi- 
cal proof in Ref. (jlll ) is not applicable directly to nuclei 
because nuclei are too dense to admit a semiclassical ap- 
proximation. 

Quantum chaos has been studied intensely during the 
last 20 years or so ()12l ). These studies have focused 
mainly on dynamical systems with few degrees of free- 
dom. The Sinai billiard is a prime example. Here chaos 



is due to the fact that the inner circle and the outer rect- 
angle possess incommensurate symmetries. In contradis- 
tinction, the atomic nucleus is a many-body system. We 
will show that in such a system, the dynamical features 
which cause chaos are distinctly different from the ones 
which cause chaos in few-degrees-of-freedom systems. 

G. Chaos in Nuclei 

Now we are in the position to return to the data shown 
in Fig. [3l Combining the information on the excited lev- 
els in 239 U obtained from these data with that of many 
similar experiments on the scattering of slow neutrons 
and of protons at the Coulomb barrier by other nuclei, 
Haq, Pandey, and Bohigas (13) formed the "nuclear data 
ensemble". It consists of 1726 spacings. These were used 
to determine the NNS distribution and the A3 statistic. 
The results are shown in Figs. [8] and [9j The solid lines la- 
beled "Poisson" correspond to totally uncorrelated levels 
and, thus, to integrable systems. The agreement with the 
GOE prediction is impressive. We conclude that at neu- 
tron threshold (and at the Coulomb barrier for protons) 
nuclei display chaotic motion. The deep physical insight 
of N. Bohr, who conceived the idea of the compound 
nucleus and who designed the wooden model of Fig. 2J 
is borne out by the data: Idealized models of billiards 
are prime examples of chaotic motion, and nuclear lev- 
els at neutron threshold display chaotic behavior. Since 
1983, there has been growing evidence for chaotic mo- 
tion in nuclei also at lower excitation energies, although 
the statistics of such data is usually not as good as near 
neutron threshold. 



to 



0.5 





1 

^ Poisson 


■ — r -- 

NOE 






1726 spacings 






GOE 









If 1 1 

12 3 



FIG. 8 The NNS distribution for the nuclear data ensemble 
(histogram) and the GOE prediction (solid line). Taken from 
Ref. (pj). 

In conclusion, there is strong evidence for chaotic mo- 
tion in nuclei. How can this evidence be reconciled with 
the extant information on nuclear dynamics (the shell 
model)? And - given the fact that GOE spectra carry 
zero information content beyond the mean level spacing 
- which is the information content of nuclear spectra in 
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FIG. 9 The A3 for the nuclear data ensemble (data points) 
and the GOE prediction (solid line). Taken from Ref. (|13h . 



the regime of chaotic motion? These are the questions 
we address in the remainder of this paper. 



titative agreement with experimental spectra. Again, we 
will confine ourselves for simplicity to two-body forces. 
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II. DYNAMICAL ASPECTS 

A. Chaos in the Nuclear Shell Model 



FIG. 10 Level sequence in the nuclear shell model, 
from Ref. dU 



Taken 



The shell model is the universal model for the struc- 
ture of atoms and nuclei. In nuclei, it consists of a cen- 
tral potential U with a strong spin-orbit interaction (the 
"mean field" ) and a fairly weak "residual interaction" V 
which accounts for the remaining part of the nucleon- 
nucleon interaction. The central potential U defines the 
major shells shown in Fig.[10j The ls-shell, the lp-shell, 
the 2sl<i-shell, etc. in spectroscopic notation, the latter 
with the subshells S1/2, ^3/2 an d ^5/2 where the indices 
denote the nucleon spin. The subshells have different 
single-particle energies and the states within a major 
shell are, therefore, not totally degenerate. Still, large 
degeneracies remain when the shell contains more than 
a single nucleon. These degeneracies are lifted by the 
residual interaction V. That interaction conserves parity 
and is rotationally invariant and invariant under time re- 
versal. It is weak in the sense that it produces configura- 
tion mixing primarily within the same major shell. More 
precisely: The magnitude of a typical matrix element is 
small compared with the spacing of the centroids of ad- 
jacent major shells but comparable with that of adjacent 
subshells in the same major shell. In all that follows we 
will, therefore, consider only a single major shell. And 
the residual interaction is predominantly a two-body in- 
teraction, V = ^2 a< pV(r a ,ri3) although there is evi- 
dence that three-body forces are needed to obtain quan- 



The Hamiltonian of the shell model is 
H = E ("|^ A « + + E V (^,r ) . (3) 

a ^ ' a</3 

Here, f a denotes the coordinates and spin/isospin degrees 
of freedom of the nucleon labeled a. The single-particle 
part of H has eigenstates \j) and eigenvalues Sj. We use 
second quantization and write H as 



H 



4 ^ 

ijkl 



v(ij]kl)a\a]aia k . (4) 



Here, Oj creates a nucleon in a single-particle state j. The 
two-body matrix elements v(ij; hi) represent the residual 
interaction. In writing Eq. (|4]) we have not paid atten- 
tion to conserved quantum numbers like spin and isospin. 
This was done in order to keep the notation simple. More 
details concerning this aspect are given in Section III.BI 

In the shell model it is assumed that all shells but one 
are completely filled, and attention is focused on that last 
shell (the "valence shell"). Distributing the valence nu- 
cleons over the various subshells in the valence shell, one 
constructs a basis of antisymmetrized many-body states 
\Jjj). We use a short-hand notation where J stands for 
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total spin, total isospin, and parity, while /i is a running 
index. These are used to calculate the matrix elements 

H^(J) = (Jfi\H\Ju) (5) 

of the Hamiltonian (|4j) in the subspace with quantum 
numbers J. Diagonalization of H^ V (J) yields the eigen- 
values and eigenfunctions. The latter are used to calcu- 
late various transition matrix elements. This scheme has 
been used with considerable success in many parts of the 
periodic table. Exceptions are mainly deformed nuclei 
which occur in the middle between large major shells. 
When we refer to Eq. (jH) in the text further below, we 
always do so with the understanding that that equation 
refers to the valence shell only. 

In the sequel we focus attention to the 2s ld-shell (in 
short: the sd-shell). The ls-shell and the lp-shell are 
filled, taking a total of 16 nucleons. Thus, the sd-shell de- 
scribes nuclei with mass numbers between 16 (0) and 40 
(Ca). We do so for practical reasons: The dimensions of 
the Hamiltonian matrices in the sd-shell are quite man- 
ageable (maximal dimensions are of the order of 10 3 ) 
while much larger numbers occur in higher shells. We 
sometimes consider also a single j-shell with half-integer 
j. Although not realistic in the framework of the nuclear 
shell model, this idealization is useful for purposes of ori- 
entation. 

Which are the fluctuation properties of the resulting 
spectra? That question was studied, for instance, in 
Ref. (15). The authors used standard single-particle en- 
ergies and a standard two-body residual interaction to 
calculate the positions of the states with spin J = 2 
and isospin T = for 12 nucleons in the sd-shell. The 
matrix dimension is 874. In Figs. [11] and [T2J the NNS 
spacing distribution and the A3 statistic for these lev- 
els are compared with GOE predictions. The agreement 
is very good. (It would be desirable to confirm this re- 
sult experimentally. Unfortunately, the data needed are 
not available: There is no more than a handful of lev- 
els known for each value of spin and isospin; such levels 
cluster in the vicinity of the ground state). In Ref. ([TBI ) 
numerous further results (on thermodynamic properties, 
on the statistics of eigenfunctions etc.) were described. 
While it was found that the input used in shell-model 
calculations (the distribution of many-body matrix el- 
ements) differs from that used in the GOE, the results 
agree to a very large extent with GOE predictions. 

The result displayed in Figs. [TT1 and fT2l and the other 
results obtained on chaos in shell-model calculations are 
surprising and call for a deeper analysis. Indeed, the 
Hamiltonian of the shell model is quite different from 
a GOE matrix. And yet, it is able to produce GOE- 
like spectral fluctuations! It does not seem unreasonable 
to speculate that the shell model would also be able to 
account for the GOE-like spectral fluctuations shown in 
Figs. [8] and [9] (this cannot be checked because of the huge 
dimensions of the matrices involved). At the same time, 
the shell model accounts well for some of the regular fea- 
tures observed experimentally like magic numbers and 
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FIG. 11 NNS spacing distribution for the J = 2, T = states 
of 12 nucleons in the sd-shell (histogram) compared to the 
GOE (line). Taken from Ref. 
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FIG. 12 A3 statistic for the J — 2, T — states of 12 nucleons 
in the sd-shell (data points) compared to the GOE (solid 
line). Taken from Ref. (0). 



spectroscopic data on low-lying states. The questions 
raised at the end of Subsection II. Gl must, therefore, be 
reformulated as follows: Are the results shown in Figs.fTTl 
and [12] a generic property of the shell model or are they 
due to very special circumstances? In other words: How 
is chaos produced in the shell model? And what is the 
information content of shell-model spectra? A further 
question would be: How does the shell model manage to 
produce both regularity and chaos? We will not address 
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that last question. Suffice it to say that regularity is 
mostly seen when we consider states with different quan- 
tum numbers like in Fig. [21 while chaos is manifest in long 
sequences of states carrying identical quantum numbers. 
The way we use them here, regularity and chaos are not 
neccessarily contradictory concepts. 



B. Two-body Random Ensemble (TBRE) 

To answer these questions, it is neccessary to go some- 
what more deeply into the details of shell-model calcula- 
tions. The single-particle states in the shell model carry 
the half-integer spin quantum number j. The two-body 
interaction scatters a pair of nucleons in states (ji, J2) 
into states (j'3, J4) which may or may not be identical to 
(ji, 32)- The two-body interaction conserves total spin, 
total isospin, and parity. Coupling the initial and the fi- 
nal pair of states to total spin s (and suppressing the anal- 
ogous coupling for isospin for simplicity) , we write the an- 
tisymmetrized reduced matrix elements of the two-body 
interaction in the form (jsj4s\V\jij2s) or, more simply, 
as v(a) where the index a labels different two-body ma- 
trix elements. (The word "different" refers to matrix ele- 
ments not connected by symmetry properties.) Counting 
(which must include spin and isospin) shows that in the 
sd-shell, a ranges from 1 to 63, and in a single j-shell 
with half-integer j, from 1 to j + 1/2. 

Chaos has to do with the complete mixing of states in 
Hilbert space. How does the residual interaction accom- 
plish such mixing? We focus attention on the two-body 
interaction and assume that all subshells belonging to a 
major shell are degenerate. (It is intuitively obvious and 
has been shown numerically that the mixing of shell- 
model states becomes weaker as the differences between 
single-particle energies of subshells become comparable 
to the magnitudes of the v(a).) We accordingly omit the 
single-particle energies in Eq. (jlj) altogether. The Hamil- 
tonian H^ v of Eq. ([5]) is then completely determined by 
the two-body interaction in Eq. (jH). Thus, it has the 
form 

H^(J) = J2Ha) C^(J,a) . (6) 

a 

The sum over a is equivalent to the sum over {ijkl} in 
Eq. (j4j). Each term in the sum in Eq. (|6j) is the prod- 
uct of two contributions. The v(a) stand for the matrix 
elements denoted by v(ij; kl) in Eq. ((4]) (except that we 
have now paid proper attention to the coupling of angu- 
lar momenta) and represent the specific features of the 
two-body interaction. Except for angular momentum 
coupling coefficients, the elements of the coefficient ma- 
trices C fIU ( J, a) are the matrix elements of the operators 
a\a^aiak in Eq. ((4]). These matrices are determined by 
the major shell we are in, by the coupling scheme we have 
chosen to construct the states | J/i), and by the specific 
two-body interaction operator labeled a which defines 



the matrix element v(a). The matrices C^ v (J,a) repre- 
sent the geometry and symmetries of the shell model and 
are generic in the sense that they are independent of the 
actual choice of the residual interaction. The seemingly 
trivial Eq. ((6]) actually contains very useful information 
about chaos in the shell model. 

To answer the question "How does the residual interac- 
tion cause chaos?" in the most general terms, we wish to 
study the Hamiltonian ((6]) for the most general two-body 
interaction. This is accomplished by assuming that the 
elements v(a) of the two-body interaction are Gaussian- 
distributed random variables with mean value zero and 
a common second moment. The value of that moment 
is irrelevant since it determines only the overall scale of 
the spectrum. With this assumption, the matrices (|6|) 
form an ensemble of Gaussian-distributed random ma- 
trices. This is the two-body random ensemble (TBRE) 
of the shell model orig inally introduced in the 1970s by 
French and Wong (|l6l ) and by Bohigas and Flores (fl7l). 
The TBRE is tailored to the shell model and, therefore, 
much more realistic than the GOE. The ensemble be- 
ing defined in terms of an integration over the Gaussian 
variables v(a), statements derived by averaging over the 
ensemble hold for all members of the ensemble with the 
exception of a set of measure zero with respect to that 
integration measure. 

Within the TBRE, the v(a) play a minor role only: 
They determine the specific linear combination of the 
matrices 6^(7, a) which forms the Hamiltonian. If the 
Hamiltonian causes chaos for almost all choices of the 
v(a), that property must be inherent in the matrices 
C^y(J^a). Therefore, an analytic theory has to address 
these matrices as the fundamental building blocks of the 
TBRE. Unfortunately, such a theory is in its infancy 
as yet. An analytically more accesible ensemble is pro- 
vided by the embedded Gaussian orthogonal ensemble 
EGOE(/c) (|18l). This is defined as a random matrix 
ensemble of m Fermions occupying / degenerate single- 
particle levels with Gaussian-distributed random /c-body 
interaction matrix elements, but assumes no further sym- 
metries like spin and isospin. For k = m, the EGOE(/c) 
becomes the GOE, and for k = 2, it becomes the analogue 
of the TBRE without quantum numbers. For recent re- 
views, we refer the reader to Refs. ([19c 1201 ). Less is kown 
for the TBRE. In the sequel, we describe some insights 
that have been gained over the last few years. 



C. Comparison GOE-TBRE 

The GOE has three important properties: It is in- 
variant under orthogonal transformations (hence, ana- 
lytically tractable), it is universal, and it is ergodic. The 
TBRE probably does not share any of these properties. 
The set of matrices C flu ( J, a) is fixed. A unitary transfor- 
mation of all matrices generates another representation 
of the ensemble but not another member of the ensem- 
ble. Therefore, the TBRE is not unitarily invariant. It 
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is not clear how a non-Gaussian distribution of the v(a) 
would affect spectral fluctuation properties of the TBRE. 
In the case of the GOE, local fluctuation properties and 
global spectral properties become separated in the limit 
TV — > oo. That separation is at the root of universality 
(local fluctuation properties do not depend on the form 
of the distribution of the matrix elements). By defini- 
tion, the TBRE is linked to a specific shell; its matrices 
have finite dimension. Similarly for ergodicity: For the 
GOE, that property is proved by showing that correla- 
tion functions vanish with increasing distance of their en- 
ergy arguments. The proof uses the limit N — > oo. This 
limit does not exist in the TBRE. One is tempted to ask: 
Why bother with the TBRE? The answer is: The TBRE 
is more realistic than the GOE. Moreover, it might be 
possible to study the TBRE analytically for the case of 
a single j-shell and to prove universality and ergodicity 
in the limit j — > oo. 

Except for the mirror symmetry about the main digo- 
nal, every element of a GOE Hamiltonian matrix stands 
for an independent random variable. For N ^> 1, the 
number N(N+ 1)/2 of such variables is much larger than 
the matrix dimension N. The number n of independent 
random variables in the TBRE is typically small in com- 
parison with the matrix dimension. For the case of a sin- 
gle j-shell, n grows linearly with j while the typical ma- 
trix dimension grows exponentially with j. In the GOE, 
the analogues of the matrices C^(J^a) exist. These are 
the TV (TV + l)/2 matrices G M which either have a unit 
element somewhere in the main diagonal and zeros ev- 
erywhere else, or have a unit element somewhere above 
the main diagonal, its mirror image below, and zeros ev- 
erywhere else. The set {G M } forms a complete basis for 
real and symmetric matrices. In contradistinction, the 
matrices C MZy (J, a) do not form such a complete set. To 
be sure, every matrix C M ^(J, a) may be thought of as a 
linear combination of the G M . But the number of matri- 
ces C fjLU ( J, a) is typically much smaller than N(N+ 1)/2. 
Therefore, many other linear combinations of the G M ex- 
ist which are linearly independent of the C^ v (J,a) and 
which do not occur in the TBRE. The TBRE may be 
negatively defined by constraining all such linear combi- 
nations to be zero. 

We are going to present evidence that the matrices 
C/xi/( J, a) are the agents which cause chaos in the TBRE. 
In anticipation of this result, we give a qualitative argu- 
ment why every matrix C^ y {J : a) may be thought of as 
a single representation of the GOE. In constructing the 
many-body basis states |J/i), we may proceed as fol- 
lows. We form Slater determinants by distributing m 
nucleons over the single-particle states of a major shell. 
These are antisymmetric by construction but do not pos- 
sess good quantum numbers J. States with good J are 
obtained as linear combinations of such Slater determi- 
nants, a typical state | J/i) consisting of a linear combi- 
nation of many determinants. Although the two-body 
interaction has non-vanishing matrix elements only be- 
tween determinants that differ in the occupation num- 



bers of not more than two single-particle states, that 
fact makes sure that the matrices C^(J^a) are densely 
filled. Moreover, the coefficients in the linear combina- 
tions are determined by angular-momentum algebra, i.e., 
contain Clebsch-Gordan coefficients, Racah coefficients, 
and coefficients of fractional parentage. Although indi- 
vidually determined group-theoretically, these quantities 
combine to make the coefficients of the linear combina- 
tions almost random numbers. This stochastic aspect of 
the shell model has been emphasized by Zelevinsky and 
collaborators who employed the term "geometric chaotic- 
ity" ©. 

Thus, the matrices C^(J^a) are dense and their ele- 
ments are close to being random. Our qualitative argu- 
ment is subject to two caveats. First, the dimension of 
the matrices C^ v (J,a) is always finite. Second, the ma- 
trices C Ml/ (J, a) are not dense everywhere and may have 
a substructure with blocks that are completely empty. 
Consider, for instance, the sd-shell and the two-body 
operator ao which scatters a pair of nucleons in states 
(^3/2^3/2) m ^o states (^5/2^5/2)- That operator has 
vanishing matrix elements between all many-body states 
|J/i) which are constructed by filling only the s 1 / 2 shell 
and the d 3 / 2 shell with m nucleons. The same statement 
applies to the matrix (7^(7, ao). 



III. PROPERTIES OF THE TBRE 

In this Section we display several properties of the 
TBRE which are either fundamentally different from 
those of the GOE or have no counterpart in the lat- 
ter. The material in this Section is largely taken from 
Refs. (|2ll : I22I ; I23I ). The numerical examples are all ob- 
tained in the sd-shell or in a single j-shell. 



A. How the matrices C^(J,a) mix the states 

The sd-shell has subshells S1/2, ^5/2? ^3/2- Let rti with 
i = 1,2,3 denote the occupation numbers of these three 
subshells, and m = J2 i rti the total number of nucleons 
in the sd-shell. The many-body basis states | Jfi) can be 
ordered in blocks with fixed values of {771,772,773}. For 
m = 12 there exist 41 such partitions. The two-body ma- 
trix elements v(a) belong to one of three classes: (i) They 
leave the partition unchanged so that {77-1,77,2,77,3} — > 
{771,772,773} (28 two-body matrix elements); (ii) they 
change the partition by moving one nucleon from one 
subshell to another while leaving the second nucleon in 
its subshell so that {77,1,77,2, 773} — > {m + 1, 772 — 1, 773} etc. 
(22 two-body matrix elements); (iii) they change the par- 
tition by moving both nucleons from one subshell to an- 
other so that {ni, 772, 773} —> {ni — 1, 772 — 1, 773 + 2} cyclic 
or to { 77 1 — 2, 772 + 1, 773 + 1} cyclic (13 two-body matrix 
elements). The same classification applies to the matri- 
ces C^ v {J^a) since each of these and the corresponding 
v(a) contain the same two-body operator. Therefore, 



10 



the C MZ ,(J, a) acquire block structure. For 12 nucleons 
coupled to J = 0, T — this is shown in Fig. [T3l where 
the 839 states are ordered according to the sizes of the 
blocks to which they belong. 




100 200 300 400 500 600 700 800 

FIG. 13 (color online) The block structure of the matrices 
C MZ /(J, a) for m = 12 nucleons in the sd-shell coupled to 
total spin J = and total isospin T = 0. Red: No change of 
partition. Green: Change of partition by moving one nucleon. 
Blue: Change of partition by moving two nucleons. White: 
There is no two-body matrix element that connects the states 
in question. Taken from Ref. 

The block structure is clearly a consequence of the shell 
model with its subshells. It is bound to occur in every 
major shell. The structure is completely absent in the 
GOE where the matrices are densely filled. It is clear 
that a single matrix C^(J, a) cannot mix the states com- 
pletely. Such mixing can be accomplished only by a linear 
combination of most or all of the C fjLU ( J, a) as in Eq. (|6|). 
Put differently, it is clear that a single or a few non- 
vanishing matrix elements v(a) cannot lead to a com- 
plete mixing of the basis states, i.e., to chaos. We expect 
that we need most or all of the v(a) for that purpose. 
This observation illustrates the caveat "with the excep- 
tion of a set of measure zero" used in Subsection III.B1 
Putting one or several of the u(a)'s identically equal to 
zero reduces the set {v(a)} to a set of measure zero with 
respect to the integration over all the variables v(a). The 
observation also sheds light on and reinforces the concept 
"geometric chaoticity" (|15f ). Indeed, the mixing is both, 
made possible and limited, by geometric constraints that 
result from the coupling of single-particle states to many- 
body states with good quantum numbers. 

A few numbers may illustrate the way in which the ma- 
trices Cy, v { J, a) fill their individual blocks. As a measure 
for such filling, we use the number of C MI/ (J, a) which 
contribute on average to a given matrix element. When 
the partition is not changed (red blocks in Figure [T3]h 



that number is 25.2 ± 4.0 (out of a total of 28 matrices). 
Within the red blocks, mixing is thus very thorough. This 
statement is supported by another measure (inverse par- 
ticipation ratios) not shown here. For the green blocks, 
the corresponding number is 7.2 ±1.4 (out of a total of 
22 matrices) and for the blue blocks, 2.0 ± 0.9 (out of 
a total of 13 matrices). Similar figures are obtained in 
the case of a single j-shell. The wealth of available nu- 
merical data (see, for instance, Figs. [11] and [12]) shows 
that there is chaos for some choices of the v(a). This is 
not an accidental feature related to a particular choice of 
the residual interaction. Indeed, the numbers just given 
show that the mixing of states within each of the diag- 
onal (red) blocks is virtually complete for almost every 
choice of the residual interaction. This fact ensures that 
a few non-diagonal matrix elements suffice to also mix 
different diagonal blocks very efficiently with each other. 
Thus, the seemingly small numbers cited for the green 
and blue blocks are actually quite sufficient to accomplish 
complete mixing of all states, even though some (white) 
blocks carry zero entries. As a result, we conclude from 
the structure of the matrices C^(J^a) that chaos is a 
generic structural property of the sd-shell. Since other 
major shells are structurally similar to the sd-shell, we 
conclude that chaos is generic in the shell model at large. 
We formulate this statement with one proviso: We have 
neglected the differences in the single-particle energies 
pertaining to different subshells. In realistic cases, our 
statement applies only if the magnitude of the matrix el- 
ements of the residual interaction is sufficiently large to 
mix different subshells sufficiently strongly. That is the 
case in practice. 



While different subshells are mixed quite thoroughly 
and chaos is expected to prevail within every major shell, 
different major shells do largely retain their identity be- 
cause the magnitudes of the v(a) are typically small com- 
pared to the spacings of major shells. This is demon- 
strated by experimental facts on the behavior of the neu- 
tron strength function or of spectroscopic factors which 
unfortunately cannot be presented and explained in the 
present framework. 



We reiterate the remark that the strong mixing within 
a major shell does not preclude the existence of individ- 
ual states at higher excitation energies that can be iden- 
tified as distinct modes of motion of the nucleus (like the 
giant-dipole resonances). Deformed nuclei display collec- 
tive motion like rotational bands and, in the ground-state 
domain, deviations from Wigner-Dyson spectral fluctu- 
ations. This last fact has been extensively discussed in 
the literature (see Ref. (24) and references to earlier work 
therein). At this point, it is not clear how such features 
can be reconciled with the view of chaos in nuclei pre- 
sented in this paper. 
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B. Information content of nuclear spectra 

GOE spectra do not carry any information content. 
Spectral fluctuations in nuclei often agree with GOE 
predictions. Are such spectra void of physical informa- 
tion? We answer this provocative question in quantita- 
tive terms (|23l ) but begin with a qualitative considera- 
tion. 

The Hamiltonian ([6]) contains the matrix elements v(a) 
as parameters and is otherwise fixed. The number of 
these parameters is typically small compared with the 
matrix dimension. Therefore, a small number of data 
points (for instance, energies of levels with fixed spin J) 
suffices, in principle, to determine H^ V (J). The Hamil- 
tonian matrices of the levels with spins J' ^ J are gov- 
erned by the same set of matrix elements v(a) (only 
the matrices C^^J'^a) differ in form and dimension) 
and are, therefore, also known once the v(a) are deter- 
mined. Hamiltonians describing other nuclei (different 
mass numbers) pertaining to the same major shell are 
likewise governed by the set {v(a)} and, thus, known too. 
(Again a caveat is needed: The set {v(a)} which empiri- 
cally describes the data best may change with mass num- 
ber). Thus, the situation is radically different from that 
of the GOE because the matrices C^ u {J^a) (the agents 
of chaos) are fixed by the shell model. A small number of 
parameters determines spectra of many J values in many 
nuclei. TBRE spectra do carry information. 

We turn to the quantitative question: How reliably can 
this information be deduced from data? To this end we 
cast the Hamiltonian (j6]) in another form. We seek a 
decomposition of the Hamiltonian into a sum of matri- 
ces that can be ordered with respect to their "relevance" 
or "importance". The decomposition (|6|) is not useful 
for this purpose: The matrices C^ v (J,a) exhibit some 
degree of linear independence but a stronger metric con- 
cept is called for. For this purpose, we introduce the 
canonical scalar product for matrices, and switch from 
the matrices C^ u (J,a) to new matrices B^ v {J^a) that 
form an orthonormal basis set. With d(J) the dimension 
of the Hilbert space spanned by the state vectors |J/i), 
we define the real and symmetric overlap matrices 

S a p(J) = d-\j) Trace[C(J,a)C(J,/3)] . (7) 

These are diagonalized by orthogonal transformations 
O(J), 

{O(J)S(J)[O(J)] T } a0 = sl(J)5 a0 (8) 

where T denotes the transpose. The matrices S a p(J) are 
positive semi-definite so that s^(J) > 0. The eigenvalues 
are arranged by decreasing magnitude, s\{J) > s|(J) > 
. . . > 0. The real roots s a (J) of the s 2 (J)'s are chosen 
positive or zero. If one or several eigenvalues vanish, we 
conclude from Eqs. JSJ) and that there exist one or 
several linear combinations of the matrices C(J, a) that 
vanish identically. For the a±(J) non-vanishing eigenval- 



ues 4(J), we define 



B^(J,a) = 




,...,ai(J) . 



(9) 

Except for possible degeneracies among the eigenvalues 
the matrices B^ u (J,a) are defined uniquely. By 
construction, these matrices are orthonormal with re- 
spect to the trace, 

d-\j) Trace[£(J,a)£(J,/3)] = 5 a(3 • (10) 

Written in terms of the matrices B (J, a), the Hamilto- 
nian takes the form 

ai 

H^(J) = Y, w (J^)s a (J)B^(J,a) . (11) 

We note that the eigenvalues s^(J) indicate the weight 
and relevance of the corresponding basis matrix B(J,a) 
in the construction of the Hamiltonian {Hj). This weight 
has purely geometric origin. The dynamical weight 
w(J,a) stems from the two-body interaction and is given 
by 

w(J, a)=^2 Oapv(p) , a = 1, 2, . . . , ai . (12) 

The w(J,a) have the same distribution as the v(a): 
They are Gaussian-distributed random variables with 
mean values zero and a common second moment. The 
form ([TT]) of the Hamiltonian with the definitions © 
and ([T2]) is completely equivalent to the original expres- 
sion (|6]) . It allows us to quantify the information content 
of nuclear spectra. If for some value of J one or several 
eigenvalues s 2 a (J) vanish, the number a\{J) of variables 
w(J,a) is smaller than the number of elements v(a). For 
that value of J it is then impossible to determine more 
than ai( J) linear combinations of the v(a) from data. We 
shall see that such vanishing eigenvalues always occur. 
Moreover, if some eigenvalues are significantly smaller 
than the leading ones, then Eq. ([TT]) shows that their 
influence on the spectrum is small (we recall that the 
matrices B(J,a) obey Eq. ([TO]) ), and it will be difficult 
to deduce their values from data. We conclude that it 
is important to determine the eigenvalues s^(J). That 
can be done without prior knowledge of the form of the 
residual interaction (i.e., of the values of the v(a)) since 
the eigenvalues are determined by the matrices C fIU ( J, a), 
i.e., by the shell model itself. 

In Refs. ([22I : I23I ), the distribution of the roots s a (J) 
of the eigenvalues was given for a single j-shell with 
j = 19/2 and for some nuclei in the sd-shell. All these 
distributions are very similar. By way of example, we 
reproduce in Fig. [14] the graph for the states with isospin 
zero in the nucleus 24 Mg. 

We note that all roots change smoothly with J. The 
largest, si(J), is significantly larger than all others. For 



12 



10 




i i i i i i i i i i_ 

5 10 



FIG. 14 Roots s a (J) of the eigenvalues s 2 a (J) denned in 
Eq. dSJ) versus total spin J for the T = states in 24 Mg. 
Taken from Ref. (jH). 

the case of the single j-shell, it can be shown analyti- 
cally that the corresponding linear combination B fIU ( J, 1) 
of the matrices C^(J^a) is approximately equal to the 
matrix representation of the monopole operator which in 
turn is approximately equal to the unit matrix. The nu- 
merical results indicate that these statements are generic; 
they hold also for 24 Mg. Hence, si(J) determines es- 
sentially the centroid of the spectrum of states with 
spin J and has little influence on spectral fluctuations. 
The latter are largely determined by the remaining non- 
vanishing eigenvalues and associated matrices 3^(3, a). 
Because of the orthogonality relation ([TO]) , the latter are 
almost traceless. 

The root si(J) is followed in magnitude by four eigen- 
values that are distinct from the remaining set. We ob- 
serve that some eigenvalues in that set vanish for large 
values of J (this is because the matrix dimension shrinks 
with increasing J and eventually becomes too small to 
support a large number of eigenvalues) and that at least 
one eigenvalue vanishes for all values of J. This is because 
the operator J 2 for total spin J is a two-body operator. 
Thus, the matrix representation of the term J 2 — J( J+ 1) 
can be written as a linear combination of the matrices 
Cnv(J,ot). This linear combination vanishes identically. 

Figure [TH shows that it is difficult to determine those 
w{ J, a)'s from data which belong to the smallest eigenval- 
ues, and that the problem cannot be solved by combining 
data for different spin values. Practitioners of the shell 
model avoid the difficulty by combining a fit to data with 
ab initio values of the v(a) obtained from many-body 
theory, see, for instance, Refs (j25l; l26l). 

C. Preponderance of ground states with spin zero 

In 1998, Johnson, Bertsch, and Dean (|27l ) reported 
that in even-even nuclei (nuclei with even proton and 



neutron numbers), the TBRE yields ground states with 
spin zero much more frequently than corresponds to the 
fraction of spin-zero states in the model space. Subse- 
quent work (|28l ; I29I ) showed that similar regularities ex- 
ist in bosonic and electronic many-body systems. The 
phenomenon has received intense attention, and we re- 
fer the reader to the reviews (|30l ; l3~fl ) . The problem has 
been understood quantitatively for bosonic models of the 
nucleus within mean- field calculations (32). For the nu- 
clear shell model, two approaches yield good agreement 
with numerical simulations. Zhao et al. (j33|) devised an 
algorithm that accurately predicts the fraction of ground 
states with a given spin. The other approach considers 
fluctuations of and correlations between the J-dependent 
spectral widths and also leads to semi- quantitative pre- 
dictions. It was put forward in Ref. (21), and is described 
in what follows. 

Key to the understanding of the phenomenon is the 
observation that the spectral widths <jj of states with 
different J values are correlated. The spectral widths 
are defined by 

aj = d-\j) Trace[# 2 (J)] . (13) 
From Eqs. (fTT]) and ([TO]) , we have 

a 2 J = jrw 2 (J,a)si(J) . (14) 

Spectral widths pertaining to different J depend upon 
the same random variables v(a) (the w(J,a) are linear 
functions of the v(a), see Eq. ([12]) ) and are, therefore, 
correlated. 

Before using this fact, we relate the spectral width aj 
to the position of the lowest (or highest) state with spin 
J, more precisely, to the distance Rj of the lowest (or 
highest) level in the spectrum from the origin (the spec- 
tral radius). (We do not distinguish the lowest from the 
highest state since the v(a) have random signs.) The 
spectral radius is expected to be a random variable. We 
ask for the probability pj with which, for a given value of 
J, Rj takes maximum value. The probability of finding 
a ground state with spin J is given by pj. 

We define the scaling factor r j by writing 

Rj = Tjcrj . (15) 

Eq. ([T5]) is useful in the present context because it turns 
out that (in contrast to <jj) rj hardly fluctuates with the 
v(a). For the case of 6 fermions in the j = 19/2 shell and 
for the states with spin J = 0, this is shown in the inset 
of Fig. [T5] Similarly constant (non-fluctuating) behavior 
of rj was found for other cases (different spin values, dif- 
ferent number of nucleons, nuclei in the sd-shell). There- 
fore, the fluctuations and correlations of the widths aj 
directly affect those of the spectral radii Rj. The linear 
fits to rj, obtained as shown in the inset, are used for the 
plot of rj versus J shown in Fig. [15] (and similarly for 
other cases). The overall monotonic decrease of rj with J 
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reflects the overall monotonic decrease of the dimensions 
d(J) with J. Indeed, average shell-model spectra are 
known to have approximately Gaussian shape. It is then 
intuitively obvious that rj grows with d(J). The decrease 
of r j with increasing J shown in the figure enhances the 
chance for states with small spins J to form the ground 
state. The same is true of the odd-even staggering of 
rj (which is due to the same reason and, everything else 
being equal, gives, for instance, preference to spin zero 
over spin 1). 




FIG. 15 (color online) The scaling factor rj of Eq. ([T5]) versus 
J for six fermions in the j = 19/2 shell (solid line). Inset: 
Linear fit to r j for the states with J = 0. The dots correspond 
to 900 realizations of the ensemble. Taken from Ref. (|2ll). 

It remains to study the probabilities that the spec- 
tral widths gj attain maximum values. For the case dis- 
played in Fig. [151 these are shown as red lines in Fig. [161 
The states with lowest and highest spins have the largest 
probabilities. Here the correlations play a role: For some 
choice of the v(a), the spectral width for some other spin 
value may take an unusually large value but so do, at the 
same time and for most cases, Go and 042 . The blue line 
in Fig. [16] gives the probability that the product tjgj 
attains maximum value. Because of the effect shown in 
Fig. [151 the highest spin value is suppressed, and spin 
zero wins. The result must be compared with the ac- 
tual probabilities (dots) that a state with spin J forms 
the ground state. The agreement, although not perfect, 
shows that the explanation accounts for the main features 
of the phenomenon. 

Similar results were obtained for eight Fermions in the 
j = 19/2 shell, for 20 Ne, and for 24 Mg. Figure EH shows 
the case of states with T = in 24 Mg. 



D. Correlations between spectra carrying different 
quantum numbers 

In the TBRE the preponderance of ground states with 
spin zero is, to a large extent, due to the correlations be- 
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FIG. 16 (color online) Six fermions in the j = 19/2 shell 
with 900 realizations of the TBRE. Solid red line: Numer- 
ical probabilities for aj to have maximum value. Dashed 
blue line: Numerical probabilities for rjaj to have maximum 
value. Dots: Numerical probabilities for state with spin J to 
be the ground state. Taken from Ref. (|2lh . 




FIG. 17 (color online) Same as Fig. [16] for the T = states 
in 24 Mg. Taken from Ref. (H3). 



tween spectral widths gj pertaining to different values of 
J. The correlations are caused by the fact that all widths 
gj depend on the same random variables v(a). Eq. ((6|) 
shows that the same statement holds for the Hamiltoni- 
ans H^ V (J) pertaining to different spin values and, more 
generally, for the Hamiltonians pertaining to nuclei be- 
longing to the same major shell. 

From the point of view of the shell model, these state- 
ments are not terribly surprising. A change of the resid- 
ual interaction causes simultaneous changes in the spec- 
tra of all nuclei belonging to the same major shell. In the 
framework of a statistical description, this fact is tanta- 
mount to the existence of spectral correlations. Still, it 
is remarkable that such correlations have never been ad- 
dressed in the framework of RMT until quite recently. It 



14 



300\ 




40 f 



-40\ 
-45 




E/v 



W/v 



45 



-70 



FIG. 18 Correlations between the level densities for states with spin and spin 2 (both with T — 0) in 24 Mg. Left panel: Mean 
value of the product of the two level densities. Middle panel: Product of the mean level densities. Right panel: The correlator 
of Eq. (fTT)) . (From 400 realizations of the ensemble). Taken from Ref. (p3l ). 



is, in fact, not obvious how to model such correlations 
in canonical RMT. The standard approach consists of 
assuming that correlations do not exist. This assump- 
tion is usually stated explicitly in the statistical theory 
of nuclear reactions. Let S a b(E, J) be an element of the 
statistical scattering matrix for scattering from channel 
a into channel b at energy E and total spin J. It is 
assumed that elements pertaining to different J values 
are uncorrelated. This assumption implies the symmetry 
of compound-nucleus scattering cross sections about 90 
degrees scattering angle in the center-of-mass system. 

It is obviously of interest to test these standard as- 
sumptions of RMT using the TBRE. While tests for cor- 
relations of the statistical scattering matrix have appar- 
ently not yet been performed, tests of spectral correla- 
tions for levels both in the same nucleus carrying different 
quantum numbers, and in different nuclei do exist (|23l). 
The level density p(E, J) for levels with spin J and en- 



J 



In Fig. [19] we show the same quantities for the J = 
0, T = states in 22 Ne and in 24 Mg. Here the correlator 



ergies E^(J) (where \i is a running index) is given by 



p(E,J) = ^5(E-E^J)) . 



(16) 



As a measure of spectral correlations, we use the corre- 
lator 



C(E u J i; E 2 ,J 2 ) = p(E 1 ,J 1 )p(E 2 ,J 2 ) 



J 



-piE^Ji) p(E 2 ,J 2 ) . (17) 

In Fig. [18] we show the two terms on the right-hand 
side of Eq. (jTTJ) and the resulting correlator versus the 
energies E\ , E 2 of the two sets of spin states for the J = 
0,T = and the J = 2,T = states in 24 Mg. The 
correlator has a maximum value of about 13 percent of 
the mean value of the product of the level densities. 



r 



has a maximum value of 6 percent. In both cases, the 
existence of spectral correlations is definitely established. 



r 



The correlators shown in Figures [18] and [19] are ob- 
tained by averaging over the ensemble and are, therefore, 
entirely theoretical. Do they have any correspondence in 
reality (where we deal with a single Hamiltonian rather 
than an ensemble)? Presently, this question cannot be 
answered using experimental data - the data set is not 
sufficient. However, shell-model calculations can serve 
as a substitute. Realistic single-particle energies and an 
optimized two-body interaction (the "Brown- Wildenthal 
interaction" (|25l)) were used to calculate the energies of 
low-lying states with spins J = and J = 2 (J = 1/2 



and J = 5/2) in a number of even-even (odd-even, re- 
spectively) nuclei in the sd-shell (|23l ). The correlations 
between nearest-neighbor level spacings of the lowest few 
states with different spins were evaluated as in Eq. ([17]) . 
the ensemble average being replaced by the running aver- 
age over the set of nuclei just mentioned. Figure [20] shows 
the result (displayed in the same form as in Figs. [TBI and 
[T9]) . Again, the existence of correlations is established. 
We observe that these are particularly pronounced for 
the lowest levels where they amount to about 10 percent. 



15 




FIG. 19 Same as Fig. 1181 except for the correlations between the level densities for states with spin and isospin in 20 Ne and 
in 24 Mg. Left panel: Mean value of the product of the two level densities. Middle panel: Product of the mean level densities. 
Right panel: The correlator of Eq. (|17[) . (From 400 realizations of the ensemble). Taken from Ref. (p3h . 




FIG. 20 Same as Fig. [18] except for the correlations between pairs of nearest-neighbor spacings of low-lying levels with different 
spins obtained by averaging over a number of nuclei in the sd-shell (see text). The indices i and j label the spacings consecutively 
starting from the lowest state. Left panel: Mean value of the product of the two spacings. Middle panel: Product of the mean 
values of the two spacings. Right panel: The correlator. Taken from Ref. d23l ). 
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Our analysis establishes the existence of correlations 
between spectra carrying different quantum numbers. All 
examples were taken from the sd-shell, however, where 
the dimensions of the Hamiltonian matrices are rather 
small (typically 10 2 to 10 3 ). The maxima of the correla- 
tors are in the 10 percent range. It is not clear how these 
results carry over to other major shells where both the 
matrix dimensions and the number of independent two- 
body matrix elements are much bigger. The statistical 
theory of nuclear reactions is mainly used in such shells. 
Therefore, definitive conclusions must wait for a further 
analysis. 



IV. SUMMARY AND CONCLUSIONS 

In the first part of this paper, we have introduced ran- 
dom matrices and the concept of (quantum) chaos. Both 
are linked by the Bohigas-Giannoni-Schmit conjecture. 
Random-matrix theory predicts fluctuation properties of 
spectra in a parameter-free fashion. In many cases, the 
fluctuation properties of nuclear spectra agree with these 
predictions: The nuclear dynamics is (partly) chaotic. 
To relate this observation to known dynamical features 
of spherical nuclei, we have, in the second part, described 
the nuclear shell model (a mean-field theory with resid- 
ual interactions). We have pointed out that shell-model 
spectra also show chaos. To explain this observation, we 
have defined the two-body random ensemble (TBRE) 
as the generic Gaussian random-matrix ensemble of the 
shell model. In part three, we have displayed a number 
of properties of that ensemble, mainly using the sd-shell. 
We have focused attention on those properties in which 
the TBRE differs from standard RMT: (i) The TBRE 
generically produces chaos. Thus, chaos in nuclei is not 
due to a particular feature of the residual interaction. 
The mechanism which mixes the basis states of Hilbert 
space differs significantly from that of canonical RMT. 
Chaos is due to the "building block" matrices C^ v {J^a). 
In each major shell, these matrices are fixed once and for 
all. They are determined entirely by intrinsic properties 
of the shell model and are independent of any partic- 
ular residual interaction, (ii) The spectra of matrices 
drawn at random from a canonical ensemble of RMT 
do not carry any information content. Because of the 
role played by the fixed matrices C MZy (J, a), this is not 
so for the TBRE. On the contrary, it is possible to ex- 
tract relevant physical information from nuclear spectra. 
We have analyzed the degree to which this can be done, 
(iii) The preponderance of ground states with spin zero 
in the TBRE is, to a large extent, due to correlations 
between widths of spectra with different spin values, (iv) 
Correlations exist more generally also between spectra 
carrying different quantum numbers and in different nu- 
clei belonging to the same major shell. The implications 
of this last statement have not yet been fully explored. 

The examples used to illustrate or prove these as- 
sertions were taken from the sd-shell or from a single 



j-shell. The line of reasoning suggests very strongly, 
however, that the spectral fluctuation properties of the 
TBRE agree with RMT predictions for all major shells. 
We have confined ourselves to a two-body residual inter- 
action and to the TBRE. It is obvious that the inclusion 
of three-body forces would enhance the tendency of the 
system towards chaotic dynamics. 

Throughout most of the paper, we have assumed that 
the single-particle energies of subshells belonging to a 
major shell are degenerate. Within this approximation, 
we have shown that quantum chaos is a generic feature of 
the nuclear shell model. This means that most choices of 
the residual interaction (but not all) will result in spectra 
showing spectral fluctuations that agree with RMT pre- 
dictions. Actually, the neglect of the differences of single- 
particle energies of subshells is not completely justified. 
The magnitude of a typical two-body matrix element 
which mixes different subshells is comparable with the 
difference of the corresponding single-particle energies. 
The resulting incomplete mixing of many-body states 
within a major shell reveals itself in slight deviations from 
RMT predictions (15). Moreover, the residual interaction 
is far too weak to thoroughly mix adjacent major shells. 
This is why shell structure remains the hallmark of nu- 
clear physics, in spite of quantum chaos. It is in that 
sense that spherical nuclei display only partial quantum 
chaos. 

Low-lying states in spherical nuclei often display reg- 
ular features. This statement is not at variance with 
our conclusions. Indeed, chaos manifests itself mainly in 
the fluctuation properties of spectra. These are defined 
in terms of statistical measures tests of which require a 
large number of levels with identical quantum numbers. 
The regular features just mentioned refer to properties 
of a small number of levels in the ground-state domain 
which carry different quantum numbers. The properties 
of these levels may be particularly sensitive to a specific 
component of the two-body interaction. Modelling the 
entire residual interaction in terms of that component 
may give rise to regular motion. 

Our discussion has been confined to spherical nuclei. 
Nuclei with mass numbers lying in the middle of large 
major shells (like the rare-earth nuclei) cannot be suc- 
cessfully described in terms of the spherical shell model. 
The dimensions of the many-body Hilbert spaces are too 
large. One uses instead collective models with a small 
number of degrees of freedom. The connection between 
these and the spherical shell model is not firmly estab- 
lished. However, quantum chaos is prevalent in these 
nuclei, too, with the exception of cases of distinct sym- 
metries (34). 

There are several open questions and directions for fu- 
ture research, (i) We are still lacking a deeper analytical 
understanding of the TBRE and its fluctuation proper- 
ties. An analytical approach must be based on proper- 
ties of the matrices denoted by C^ v {J^a) in this paper. 
While a theoretical description for shells with several sub- 
shells is probably very difficult, focus on a single j-shell 
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might simplify the problem, (ii) The TBRE predicts cor- 
relations between spectra with different quantum num- 
bers (e.g., different masses, spins, or isospins) for nuclei 
within a major shell. The expermimental verification is 
difficult due to limitations in the length and complete- 
ness of observed nuclear spectra, but other Fermi systems 
might be more accessible, (iii) The correlations between 
spectra with different quantum numbers might also affect 
the scattering matrix, more precisely: Such correlations 
might induce correlations among ^-matrix elements car- 
rying different total spin quantum numbers. The present 
analysis of fluctuating cross sections in compound nuclei 
neglects any such correlations. A better understanding 
of this problem would be highly desirable. 

Nuclear spectroscopy is a mature field with a history 
of more than 50 years. In spite of this fact it continues to 
offer great challenges. We have addressed one of them: 
The way chaos is induced by the two-body interaction of 
the shell model, an interaction which is at the same time 
responsible for the many regular features seen in nuclei. 
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